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1. Introduction 

Binary systems made of compact objects (neutron stars or black holes) that spiral 
in toward coalescence because of gravitational-radiation damping are among the most 
promising candidate sources for interferometric gravitational-wave (GW) detectors, such 
as the Laser Inteferometric Gravitational Wave Observatory (LIGO), VIRGO, GEO and 
TAMA The inspiral waveform will enter the detector frequency band during the 
last few minutes or seconds of evolution of the binary and the GW community plans to 
track the signal phase and build up the signal-to-noise ratio by integrating the signal for 
the time during which it stays in the detector bandwidth. This is achieved by filtering 
the detector output with a template which is an (approximate) theoretical copy of the 
observed signal. 

Einstein theory predicts that the radiative transverse traceless (TT) gravitational 
field hj^^ , far away from the source, is related to the motion and the structure of 
the source, at lowest order in the post-Newtonian (PN) expansion, by the quadrupole 
formula 



K ^) = ^ V^O^^raW ^Qfc„ y--)^ 

where Qij {i,j = 1,2,3) is the tracefree quadrupole moment of the source; D is the 
distance from the source; N = X /D is the unit vector from the source to the observer; 
'Pijkmi.N) is the TT projection operator onto the plane orthogonal to AT"; G is the 
Newton constant; and c is the speed of light. From Eq. (0) (and its extensions at higher 
PN orders) we see that the more precisely we know the two-body motion, the more 
accurately the PN template Hf^ will describe the "real" gravitational waveform. 

In Fig. |l] we show a typical gravitational waveform. The part of the waveform drawn 
with a continuous line is emitted during the inspiral phase when the two black holes 
are largely separated, e.g., r > lOGM/c^ where we denoted by r the radial separation 
and by M the total mass of the binary system. During the inspiral, the two black holes 
follow an adiabatic sequence of quasi-circular orbits. The equation of motion in the 
center-of-mass frame can be written schematically as 

^ = [l+0(e)+0(e2)+0(e^/2)+0(e3) + - ■ •] x \l+0{v)+0{v') + . ■ •] , (2) 

where x denotes the separation vector between the two bodies and r = |a3|. Equation 
(0) is characterized by a double expansion: in the PN parameter e ~ f^/c^ ~ GM/ (c^r), 
and in the parameter v = nii m^jM'^ ^ where mi and m2 are the masses of the two black 
holes. The parameter v ranges between (test-mass limit) and 1/4 (equal- mass case). 

It is well known that the PN expansion converges badly: as the two bodies draw 
closer, and enter the nonlinear, strong curvature phase, the motion becomes relativistic, 
e.g., w/c ~ 0.3, and it is more and more difficult to extract reliable information from 
the PN series. More specifically, when the distance between the inspiraling black holes 
shrinks to r < lOGM/c^ the PN expansion can no longer be trusted . The dashed line 
in Fig. m depicts the (less known) part of the waveform emitted during the final phase 
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Figure 1. Typical gravitational waveform emitted throughout the inspiral, plunge 
and ring-down phases. 

of evolution, when nonlinearities and strong curvature effects become important and 
nonperturbative analytical and/or numerical techniques should be used to describe it. 
This final phase includes the transition from adiabatic inspiral to plunge, beyond which 
the two-body motion is driven (almost) only by the conservative part of the dynamics. 
For nonspinning binary black holes, the plunge starts at the innermost stable circular 
orbit (ISCO) of the binary black holes. Beyond the plunge the two black holes merge, 
forming a Kerr black hole. As the system reaches the stationary Kerr state, the nonlinear 
dynamics of the merger resemble more and more the oscillations of the black-hole quasi- 
normal modes |0] and the gravitational signal will be a superposition of exponentially 
damped sinusoids (ring-down waveform). However, if the black holes carry a big spin, 
it is likely that the plunge stage will never be reached and the late dynamical phase will 
be much more complicated in this case. 

Likely, the first detection of gravitational waves with LIGO/ VIRGO interferometers 
will come from binary systems made of black holes of comparable masses, say with a 
total mass M ~ lOM© BOM©. In fact, if we assume binary black holes are originated 
by "population synthesis" f^, the estimated detection rate per year is <4 x 10^^ — 0.6 
at 100 Mpc 1^, while if globular clusters are considered as "machines" for making binary 
black holes the detection rate per year is ~ 0.04 — 0.6 at 100 Mpc 0. These numbers, 
especially the latter, are more optimistic than the estimated detection rate per year for 
neutron star binaries, which is <3 x 10~^ — 0.3 at 20 Mpc or for neutron star/black 
hole binaries, which is <4 x 10""^ - 0.6 at 43 Mpc [g. 

Although the number of cycles in the LIGO/VIRGO frequency band of gravitational 
waves emitted by comparable-mass binary black holes is not very high, on the order of 
10 — 200, however these particular sources demand a more careful analysis because the 
gravitational waves the detectors will be sensitive to are emitted during the final stages 
of inspiral where PN expansion fails. For example, in the nonspinning case the GW 
frequency at the ISCO (evaluated using the Schwarzschild ISCO) is /gw° — 220 Hz for 
M = 2QMq and ~ 167 Hz for M = BOM©, well inside the LIGO/VIRGO band. 

Moreover, black holes could carry big spins which could affect the waveforms, so if data 
analysis will be done with spinless templates there is considerable chance to miss the 
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gravitational-wave signal. 

In the next section we shall summarize what has been done in the literature to cope 
with those problems and which issues remained to be solved. 

2. Analytic methods to predict gravitational waveforms 

For spinless black hole binaries, we consider the so-called restricted waveform: h{t) = 
v'^ cos((y9Gw(^))) where v^gw = 2(^9 with ip the orbital phase and v is the invariantly 
defined velocity v = (M LpY^^ = (vtM/gw)^^^- This waveform is obtained disregarding 
all the multipolar components appearing in the gravitational waveform except the 
quadrupolar one [see Eq. (|])]. To determine in the adiabatic limit the evolution of 
the GW phase, v'gw(^)) it is sufficient to use the energy-balance equation 

relating the orbital energy function S (center-of-mass energy that is conserved in absence 
of radiation reaction) to the gravitational-ffux (or luminosity) function JF, which are 
known for quasi-circular orbits as a PN expansion in v. It is easily shown that Eq. (j^) 
is equivalent to the following system of differential equations [see, e.g., [PTI 

dt ~ M dt~ MdS{v)/dv ' ^ ^ 

whose solution provides the phasing ipcwit) during the inspiral. 

2.1. "Genuine" post-Newtonian calculations 

Let us summarize what we know about the two crucial ingredients S and JF entering 
Eq. (^). The equations of motion of two compact bodies at 2.5 PN approximation 
were ffrst derived in Refs. [0. The 3PN equations of motion have been obtained by 
two separate groups: Damour, Jaranowski and Schafer p| used the Arnowitt-Deser- 



Misner (ADM) canonical approach while Blanchet, Faye, and de Andrade [10] worked 
with the PN iteration of Einstein equations in harmonic gauge. Recently, Damour 
et al. [^ working in the ADM formalism and applying dimensional regularization, 
uniquely determined the so-called "static" parameter which entered the 3PN equations 
of motion |P, |lOl and up till then was unknown. Thus, at present time the energy 



function S is known up to 3PN order. 

The gravitational ffux emitted by compact binaries was ffrst computed at IPN order 
in Ref. |Tl|. Subsequently, it was determined at 2PN order using a formalism based on 
multipolar and post-Minkowskian approximations and independently by using a direct 
integration of the relaxed Einstein equations [0. Non-linear effects of tails at 2.5PN 
and 3.5PN orders were computed in Refs. ||T3|. More recently, the gravitational-ffux 
function for quasi-circular orbits has been derived up to 3.5PN order WM. However, at 



3PN order the gravitational-ffux function depends on an arbitrary parameter which 



is not ffxed by the regularization scheme used by the authors. 
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Although the knowledge of higher order PN corrections to S and ^ is a necessary 
ingredient to extract the phase of gravitational signals emitted by neutron star or black- 
hole/neutron-star binaries, however it is not by itself sufficient for computing GWs 
of comparable-mass binaries. As underlined above, this is due to the fact that for 
these sources LIGO/ VIRGO will detect gravitational signals emitted when the motion 
is relativistic and (genuine) perturbative PN calculations can no longer be trusted. 



2.2. Post- Newtonian resummation methods 

Certainly, the best way of extracting the GW signal emitted by comparable-mass 
binaries during the last stages of inspiral, would be to solve numerically the Einstein 
equations of a binary black hole system. Unfortunately, despite the interesting progress 
made by the numerical relativity community during the recent years JTS], |17|, |T8|, p!9|] . 



an estimate of the waveform emitted by black-hole binary has not yet been provided. 
Preliminary results for the plunge, merger and ring-down waveforms were only recently 
obtained and they use initial conditions at the ISCO which differ from PN 
predictions. To overcome this gap and tackle the delicate issue of the late dynamical 
evolution, various nonperturbative analytical approaches have been proposed ||2l|, 
24| to study the motion of two spinning and nonspinning bodies in general relativity. 

The main features of the various PN resummation methods can be summarized 
as follows: (i) they provide an analytic (gauge- invariant) resummation of the orbital 
energy function S and gravitational flux function JF, (ii) they can describe the motion 
(and provide the gravitational waveform) beyond the adiabatic approximation and (iii) 
they can in principle be extendible to higher PN orders. More importantly, they can 
be used to provide initial data for black holes just starting the plunge motion which 
can be used by the numerical relativity community to evolve the full Einstein equations 
during the merger phase. However, the resummation methods are also based on some 
assumptions hard to prove rigorously - for example, in deriving the orbital energy and 
the gravitational flux functions in the comparable-mass case, it is assumed that they 
are smooth deformations of the analogous quantities in the test-mass limit. Moreover, 
in absence of an exact solution or of experimental data we can test the robustness and 
reliability of those resummation approaches only using internal convergence tests. 

As underhned at the beginning of Sec. 0, in absence of spins, the two crucial 
ingredients necessary to extract the GW phasing are the orbital energy function S 



and the gravitational flux function J-'. In Sec. |2.2.1| we shall discuss a resummation 
approach which provides a better behaved flux-type function [^, while in Sec. |2.2.2 
we shall summarize another resummation method which provides an improved energy- 
type function ||22|, |23]. By combining these two resummation methods, Buonanno and 



Damour pSf proposed a way of describing the black hole motion beyond the adiabatic 



limit, including the transition from inspiral to plunge, whose main features are discussed 



in Sec. 2.3 
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2.2.1. Fade approximants Starting from the PN expansions oi S{v) and J-'{v) Damour, 



Iyer and Sathyaprakash |2T| proposed a new class of waveforms based on the systematic 
use of Pade resummation, || which is a standard mathematical technique used to 
accelerate the convergence of poorly converging or even divergent power series. For 
lack of space we shall discuss briefly only the Pade approximant to the flux function J-'. 

In the test mass limit (z/ — »• 0) the flux function has a simple pole at the light ring 
(f ^ = 1/3). Damour et al. |^ argued that the origin of this pole is quite general and that 
we should expect a pole singularity also in the equal-mass case. Thus, after factoring 
out this pole (defined by what they call the better behaved energy-function e{v), see 
Ref. 1^), they introduced the better behaved flux-function /(t>) = (1 — v/vpo\c) ^{v). 



Using the 2.5PN expansion of J-'{v) and applying the (near diagonal) Pade resummation 
H to f{v), they derived the f"-Pade approximant to the gravitational flux function: 

To test the reliability of the result, Damour et al. showed that in the test mass limit 
the (near diagonal) Pade approximant to J^{v; 1^ = 0), given by Eq. (^, exhibits a very 
good convergence toward the exact result, which is numerically known when u = 0. 
[See Fig. 3 in Ref. Thus, arguing that the equal- mass case can be obtained as a 



smooth deformation of the test mass limit, with deformation parameter z/, they propose 
J-p{v; v) as the best estimation of the GW flux for comparable-mass binaries. 

2.2.2. Effective- one-body reduction The resummation technique discussed in this 
section, the so-called effective-one-body (FOB) approach [^, was originally inspired 
by a similar approach introduced by Brezin, Itzykson and Zinn- Justin to study 
electromagnetically interacting two bodies. The basic idea, illustrated in Fig. 0, is 
to map the real conservative two-body dynamics up to 2PN order (see below for the 
extension at 3PN order) onto an effective one-body problem, where a test particle of 
mass /i = mim2/M, with mi, m2 the black- hole masses and M = mi + m2, moves in 
some effective background metric g'^^. This mapping has been worked out within the 
Hamilton- Jacobi formalism, by imposing that whereas the action variables of the real 
and effective description coincide, i.e. Lreai = -Z^es, Ireax = ^es, where L denotes the 
total angular momentum, and X the radial action variable, the energy axis is allowed 
to change, i^reai = f{Ecs), where / is a generic function. By applying the above rules 

J Let us assume that we know the function g{v) only through its Taylor approximant GAr(w) = 
5o + ffi + • • • + ffAf v'^ = TN[g{v)]. The idea of Pade summation ||2^ is to replace the power series 
Gjv(w) by the sequence of rational functions 

Mr / M _ ^m{v) _ Ejio«J«^ 



with M + K ^ N &Yid TM+KlPMlgiv)]] = Gn{v) = !)• For M,K ^ +oo, P*^(v) g{v). The 
choice K — M oi K = M + 1 defined the near diagonal Pade approximants. 
§ See previous footnote. 
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Figure 2. How the BOB approach matches the real two-body problem (on the left) 
and the effective one-body problem (on the right). [We defined N — I -\- L, where I is 
the radial action variable and L the angular momentum.] 

defining the mapping, it was found tliat as long as radiation-reaction effects are not taken 
into account, the effective metric is just a deformation of the Schwarzschild metric, with 
deformation parameter u = fi/M. The effective metric reads ||22|| : 



-AiR) dt" + dR^ + R^ d^^ 



(6) 



More importantly the reduction to the one-body dynamics provides the improved real 
Hamiltonian 



-T/improved 
''-real 



\ 



l + 2u 



' ni^{R,PR,P^)-lic'' 



(7) 



where {r,Pr,Pip) are the ADM coordinates used in the real description. The effective 
coordinates (i?, Pr, P^) can be related to 



transformation given at 2PN order in Ref. 
effective Hamiltonian reads: 



r,pr,p^) by a (generalized) canonical 
2l and at 3PN order in Ref. il. The 



nos{iy,R,PR,P^) 



jJ,C 



\ 



AiR) 1 + 



A{R) Pi 



+ 



p2 



(8) 



l2^c^D{r) ' ii'^c^R'^ 

We refer to the real Hamiltonian (0) as improved real Hamiltonian because being a 
(partial) PN resummation of the original badly convergent real Hamiltonian Tireai = 
'^^New + 'Hip^/c^ + 'H2vn/c^ ■ ■ ■, it should better capture the nonperturbative effects of 
the final stage of the black-hole motion. Remarkably, the mapping between the real and 
the effective Hamiltonians in Eq. (0) coincides with the mapping obtained in Ref. 



in the context of quantum electrodynamics, where these authors mapped the one-body 
relativistic Balmer formula onto the two-body energy formula. 

The EOB approach was then extended at 3PN order in Ref. |Q . The authors found 
that starting from the 3PN level there are more equations to satisfy than the number of 
free parameters appearing in the energy-map and in the effective metric. Hence, they 
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Figure 3. Comparison between the PN-expanded (on the left) and the EOB (on the 
right) binding energies for equal-mass binary. 

suggested the following two possibilities. At the price of modifying the coefficients of the 
effective metric at IPN and 2PN levels, and the energy-map (|^) as well, it is still possible 
at 3PN order to (uniquely) map the real two-body dynamics onto the dynamics of a 



test mass moving on a geodesic [see for detail Appendix A of Ref. [0]]. However, this 
solution looks quite complicated and, more importantly, it does not look very natural 
to wait to know the 3PN Hamiltonian to derive the matching at IPN and 2PN levels. 



The authors then suggested to abandon the hypothesis (used at 2PN order that 
the effective test mass moves along a geodesic and introduced in the Hamilton- Jacobi 
equation (arbitrary) higher derivatives terms which provide enough coefficients to obtain 
the matching. Because of these terms the effective 3PN Hamiltonian is not uniquely 
fixed by the matching rules defined above; the general expression is given by Eq. (3.12) 



in Ref. It is interesting to note that also in the non-geodesic case the relation 

between the effective and real Hamiltonians is still given by Eq. (|^). 

To have a qualitative understanding of the way the EOB can accelerate the PN 
convergence, we compare in Figs. ^ the PN-expanded binding energy (on the left panel) 
to the EOB-resummed binding energy (on the right panel) versus the radial separation, 
for circular orbits, in the equal-mass case, for two choices of the angular momentum L 
and at different PN orders. [To obtain the plots we used the 3PN-expanded and EOB- 



resummed Hamiltonians given by Eqs. (3.12) and (4.27) of Ref. ||2^ with ci;^ = and 



zi = = Z2. We also used the (generalized) canonical transformation between the ADM 
and the "effective" coordinates given by Eqs. (4.29)-(4.36) of Ref. |23|.] 



Note that the PN-expanded energy oscillates at the various PN orders while the 
EOB energy has a monotonic behaviour. The fractional difference between the PN- 
expanded and EOB binding energies is, e.g., for L = 3.2/iGM, 2.3% (2PN order) and 
0.24% (3PN order) at R = QGM/c^ and 17.6% (2PN order) and 2.8% (3PN oder) at 
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Figure 4. Summary of the binding energy at the ISCO predicted by post-Newtonian 
and numerical relativity calculations for nonspinning binary black holes of total mass 
15M0 -I- 15M0. 



R 



An interesting nonperturbative information of binary black hole systems is the 
presence of the ISCO defined as the solution of the equations dT-i/dr = = 9^7Y/9r^. 
In the test mass limit (Schwarzschild metric) the ISCO exists and the nonrelativistic 
energy associated to it is £f^^ = — 0.01430Mc^. If we consider the PN-expanded real 



Hamiltonian [Eq. (4.27) of Ref. ||2^] in the test mass limit and look for the ISCO, we 
find that at IPN order S™'""'^ = -O.OOTTSMc^, at 2PN order the ISCO does not exist, 
while at 3PN order ^j™"*^"^^ = — 0.01129Mc^. On the other hand, because by definition 
the EOB converges in the test mass limit to the Schwarzschild case, it automatically 
provides the correct ISCO. Thus, in the same way as Pade approximants, the EOB 
resummation method, provides by construction the right prediction in the test mass 
limit. Although we cannot prove that in the comparable-mass case the EOB approach 
is converging to the right limit, however, we can certainly say that it shows reasonable 
stability in the predictions at IPN, 2PN and 3PN orders. 

In Fig. I we summarize the binding energy at the ISCO predicted by various post- 

[51 



Newtonian ^ and numerical relativity calculations |T^, ^7\, . The equal- 

mass ISCO binding energies predicted by the EOB approach at various PN orders 
are: £^g^ = -0.01440Mc2 (IPN order), Sf^g^ = -0.01498Mc2 (2PN order) and 
£:.EOB ^ -0.01670Mc2 (3PN order). Note that the fractional difference from IPN 
to 2PN order is 4% and it increases to 10% from 2PN to 3PN order. Figure ^ shows 
that except the very recent result of Ref. |jl8|, which is much closer (probably due to 
their definition of the binary orbital frequency) to the PN estimates than the numerical 
relativity predictions, the PN and numerical relativity results |T5| , |T^, ^ differ very 
much. 
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2.3. Beyond the adiabatic approximation: transition inspiral to plunge 



Let us now introduce radiation-reaction effects. As discussed in Sec. U.2.1V Pade 



approximants |^ give a good estimate of the energy-loss rate along circular orbits 



up to 2.5 PN order. In Ref. this resummation method was combined with the 



EOB approach, and the authors deduced a system of ordinary differential equations 
which describe the late dynamical evolution of a binary-black-hole system. In spherical 
coordinates {ip, R, P^, Pr), their relevant equations are p3| : 



dR 




dPR 


Itt 


dPR ' 


dt 


dip 


rj-T/impr 


dP^ 


It 


dP^ ' 


dt 



^ dR 



Differently from Eqs. (@), the above equations go beyond the adiabatic approximation (at 
least for what concerns the conservative part of the dynamics) and can be analytically 
or numerically solved to study the transition between the adiabatic inspiral and the 
plunge. 

Let us discuss briefly the main features of this transition in the two extreme limits 



z/ ^ 1 and u = 1/A [see Ref. [^] for more details]. The case z/ -C 1 refers to binary- 
black-hole systems in which a very small black hole spirals around a supermassive black 
hole. They are typical GW sources for the Laser Interferometer Space Antenna (LISA). 
In this case, it was found ^ that the transition from adiabatic inspiral to plunge 
is sharply localized around the ISCO. Ori and Thorne ||28[ pointed out that likely LISA 



could observe the transition from inspiral to plunge with a signal-to-noise ratio of a few. 

For equal- mass binaries {u = 1/4), the radiation damping effects are important in 
an extended region on the order of A{Rc^/GM) ~ 1 above the naive (Schwarzschild) 
ISCO R = GGM/c^. The transition from inspiral to plunge is rather blurred and the 
dephasing between the full and the adiabatic waveform becomes visible somewhat before 
the naive ISCO. The plunge part of the exact waveform looks like a continuation of the 
inspiral part because the orbital motion remains quasi-circular throughout the plunge. 

Figure ^ shows the plunge, merger and ring-down part of the waveform [^] for 
equal-mass binary. The ringdown waveform contains only the mode that is damped 
more slowly, Z = 2, m = 2 at frequency Uqnm ~ 1880 (IOMq/Mbh) Hz, where Mbh 
is the mass of the final black hole formed. The dimensionless rotation parameter is 
obh = <^BH/(GMgjj) = 0.795, where we denoted the angular momentum of the final 
Kerr black hole by Jbh- The energy emitted during the plunge is ~ 0.7% of M, with a 



comparable energy loss ~ 0.7% of M during the ring-down phase |30|. This gives a total 
energy released of ~ 1.4% of M to be contrasted with the much larger value 4 — 5% of 
M recently estimated in Ref. . 



Recently Damour, Iyer and Sathyaprakash ||2^ investigated the consequences of the 
EOB waveform for LIGO/VIRGO data analysis. They derived that the GW radiation 
coming from the plunge and merger can significantly enhance the signal-to-noise ratio 
for binaries of total mass M > 30Mq. They found that the signal-to- noise ratio reaches 
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Figure 5. Plunge and ring-down gravitational waveform obtained from the EOB 
approach for equal-mass binary. 



the maximum value of ~ 8 for M ~ 80Mq at 100 Mpc. Previous estimations using 
maximally spinning binaries found that the merger dominates on inspiral. For example 
Flanagan and Hughes predicted that the energy released during the plunge is 
~ 10%M while during the ring-down phase ~ 3%M and the signal-to-noise ratio reaches 
the maximum value of ~ 40 for M ~ 2OOM0 at 100 Mpc. 



3. Open issues 

3.1. Spinning binary black holes 

The theoretical prediction of GWs from comparable-mass binaries is not only affected 
by the failure of the PN-expansion, but also by spin effects. Various studies Q 
estimated that if the binary's holes carry a spin, the evolution in time of the GW phase 
will be significantly affected by it - for example the spin can introduce modulations 
and irregularities in the gravitational waveforms. These features can become quite 
"dramatic" as long as the two spins are big and not aligned or antialigned with the 
orbital angular momentum. 

To give a rough idea of the effects we are alluding to and that we are currently 
investigating we sketch in Fig. |]the Fourier- domain phase 0(/) of the GW signal 
versus frequency for various approaches, including possible modulations due to spin 
effects. The various methods give the same prediction for the Fourier- domain phase up to 
the frequency where the PN series fails - for example for M = 20Mq, /pN-faiiure ~ 50Hz, 
while for M = SOM©, /pN-faiiurc ~ 33Hz [§, which are well inside the LIGO/VIRGO 
band. 



More recently. Levin |3J] claimed that because of spin effects the two-body dynamics 
could be affected by chaos, or more in general the dynamics, and as a consequence 
the gravitational waveform, could depend strongly on the initial conditions. However, 



Schnittman and Rasio reanalysed this issue in Ref. By calculating the divergence 
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Failure ofPN 



Mf 

Figure 6. Expected Fourier-domain phase 0(/) of the GW signal versus frequency 
for various approaches modehng the two-body dynamics. 

of nearby trajectories for a broad sample of initial conditions they concluded that the 
divergent time is much greater that the inspiral time. So even if chaos were present it 
should not affect the detection of inspiral waveforms with LIGO/VIRGO. 

To tackle these delicate issues it would be desirable to extend the resummation 
methods discussed above to the case of spinning binaries. To this respect we notice 
that the Fade resummation method has been recently extended to the gravitational 
flux function including spin-orbit and spin-spin effects Moreover, on the line of 

the EOB appoach, Damour has recently mapped the conservative dynamics of two 
spinning black holes into the one of a test particle moving in an external //-deformed 
Kerr metric. It will be very important to complete Damour's analysis by including 
radiation reaction effects and describe the late nonadiabatic dynamics of spinning black 
hole binaries. 



3.2. Possible strategy for not missing the GW signal from comparable-mass binary 
black holes 



In sections |2.2.1| and |2.2.2| we discussed the best theoretical gravitational waveforms 
currently available. These waveforms should certainly be used in detecting and 
extracting physical parameters (masses and spins) for GWs emitted by neutron-star 
binaries and neutron-star/ black hole binaries. On the other hand the direct application 
of these templates to detect GWs from comparable-mass black-hole binaries is less 
straightforward. Indeed, the approaches analysed in sections |2.2.1| and |2.2.2| are 



inevitably affected by approximations and assumptions, introduced to overtake the 
fact that the two-body problem is known only up to a certain PN order. It would 
be important for detection purposes to quantify such uncertainties and take them into 
account when building the best theoretical GW templates. If not, the risk is to miss 
the GW signal. A possibility is to deform and expand (e.g., by introducing new 
parameters), the template family generated by the resummation method having better 
PN convergent properties, in such a way that the new (higher-dimensional) template 
space could (i) include, or at least be not very far, from all the other approaches 
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characterized by worse PN convergent properties and (ii) describe signals (hopefully 
the real GW signal!) whose functional form cannot be described by any of the original 
template families. To reach this goal the ideas underlying the Fast Chirp Transform 
technique, recently proposed in Ref. EOi, deserves attention. Moreover, for comparable- 



mass binaries the inclusion of spin effects and the enlargement of the template space, 
hopefully, should not make the total number of templates huge, because the number of 
cycles expected in the LIGO/VIRGO band for nonspinning binaries is already rather 
small ~ 20—100. [The situation would be different for neutron star or neutron-star /black 
hole binaries for which in the nospinning case we expect a number of cycles on the order 
of ~ 10^ - lOM 



The new template space could be used for on-line search. When an acceptable 
signal-to-noise ratio is found, say >8, then the output signal could be re-analysed 
by templates provided, e.g., by the PN Taylor-expanded, Pade approximant and EOB 
methods, to (i) determine which approach better describe the real signal and (ii) extract 
the binary's parameters, as the masses and the spins. 
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Note added. While this manuscript was being refereed, an interesting paper of 
Blanchet appeared on the web. He compares the minimum energy of circular orbits 



for equal-mass-binary, at 3PN order, with the numerical-relativity results of Ref. 
finding good agreement. 
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